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^ '■ Abstract 

We present a theoretical and experimental study of the Fiske resonances in the current-voltage 
characteristics of "small" Josephson junctions with randomly distributed misaligned Abrikosov 
vortices. We obtained that in the presence of Abrikosov vortices the resonant interaction of elec- 
tromagnetic waves, excited inside a junction, with the ac Josephson current manifests itself by 
Fiske steps in a current-voltage characteristics even in the absence of external magnetic field. We 
found that the voltage positions of the Fiske steps are determined by a junction size , but the Fiske 
step magnitudes depend both on the density of trapped Abrikosov vortices and on their misalign- 
ment parameter. We measured the magnetic field dependence of both the amplitude of the first 
Fiske step and the Josephson critical current of low-dissipative small Nb based Josephson tunnel 
£Nj ■ junctions with artificially introduced Abrikosov vortices. A strong decay of the Josephson critical 

o , 

current and a weak non-monotonic decrease of the first Fiske step amplitude on the Abrikosov 
vortex density were observed. The experimentally observed dependencies are well described by the 

• 5-h , 

^ ■ developed theory. 

.«?; 
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I. INTRODUCTION 



It is a well known that the current-voltage (/ — V) characteristics of Josephson tunnel 
junctions can display sharp current resonances. These interesting features, called Fiske 



osephson current and standing 
l|. A long time ago the Fiske 



steps, appear as a result of resonant interaction between ac 
electromagnetic waves that can be excited in the junction 
steps have been experimentally observed in planar Josephson tunnel junctions [2j. In this 
time the complete theoretical description of the Fiske steps in the case of a junction with low 
quality factor has been done by Kulik in jjj. After that the Fiske steps have been observed in 



numerous Josephson coupled systems, e.g. Josephson junction arrays, stacked artificial 



Josephson junctions |5j, intrinsic high-T c Bi2Sr 2 CaCu20 8 +x stacked superconducting tunnel 
junctions jfj], high-T c Y Ba^CuzO-j-^ bicrystal Josephson junction [7]. Recently, the Fiske 
steps have been investigated theoretically in — tx Josephson junctions {§1 and they have 
been experimentally observed in — tt Josephson tunnel junctions with ferromagnetic barrier 

2]. 

In "small" Josephson junctions in which a junction size is less than the Josephson pene- 
tration depth Xj, the spectrum of the electromagnetic waves has a form: u(k) = Cok, where 
Co is the Swihart velocity, and the wave vectors k n are determined by a junctions size W as 

k n = nn/W, and n = 1, 2 The voltage positions of n-th Fiske steps V n are determined 

by the spectrum on(k n ) of the electromagnetic waves as []| 

v n = (i) 

2e V ' 

The magnitude of Fiske steps depends strongly on the spatial dependence of the Josephson 
phase difference between the electrodes <p{p), where p = {x,y} are the coordinates in the 
junction plane. 

In the absence of inhomogeneities a simplest way to create a coordinate dependence 
of the Josephson phase difference ip is to apply an external magnetic field parallel to the 
junction plane. In this case the Josephson critical current is strongly diminished but the 
magnitude of the Fiske step increases substantially. The magnetic field dependence of Fiske 
steps has been investigated in Josephson tunnel junctions with different shapes, namely, 
square shaped junction [loj], circular shaped junction [ll], quatric shaped junction jl2]. In 



13( the Fiske steps have been experimentally and theoretically investigated in the case of 
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the annular shaped Josephson junction both in an external magnetic field and in a magnetic 
field generated by the injection current passed along one of the annular junction electrode. 
However, even in the absence of externally applied magnetic field the Fiske resonances 
can appear in nonuniform Josephson junctions. Indeed, as was observed in Ref. jjjj], the 
randomly distributed barrier inhomogeneities (the structural fluctuations) provide necessary 
conditions for generation of Fiske steps without external magnetic field. 

From the experimental point of view the nonuniform Josephson junctions can be formed 
by inclu ding artificial inhomogeneities in the tunnel barrier 15| or by introducing Abrikosov 



vortices 



16 



17]] . Inhomogeneities of the first type suppress locally the Josephson tunneling, 



and, therefore, the spatial variation of Josephson critical current density occurs. These 
inhomogeneities may be included in the junction barrier during the fabrication process. So 
a variation of some property of the inhomogeneities (for example, their density) implies 
necessity to fabricate several Josephson junctions in which this property varies. Obviously, 
the practical realization of this kind of experiment is a rather complicated task. 

On the other hand the Abrikosov vortices can be easily trapped in the junction by ap- 
plying an external magnetic field perpendicular to the junction plane. Moreover, if pinning 
centers are randomly distributed in superconducting electrodes, the Abrikosov vortices can 
be trapped in the particular form of a misaligned vortex such that the magnetic field of the 
vortex enters and leaves the superconducting electrodes of a junction in different points (see 
Fig. 1). We notice here that an extreme case of misaligned vortices, so called "pancakes" 
vortices, naturally appears in high-T c superconductors 18[. The misaligned Abrikosov vor- 
tices can be considered as local magnetic inhomogeneities leading to an additional spatial 
variation of the Josephson phase difference f(p), and, hence, the variation of the Josephson 
critical current density over a junction area jla, O, Q, 13 . By means of field cooling process 
one can vary the density of Abrikosov vortices in the same junction. The possibility to trap 
Abrikosov vortices in a relatively easy way offers many opportunities for experimental stud- 
ies of Josephson junctions with inhomogeneities. Thus, the dependence of Josephson critical 
current on the Abrikosov vortices density and an externally applied magnetic field have been 



intensively studied 



13, 



19H22j many years ago. Recently, in Refs. [23|, |24j a Josephson 



junction has been used as a tool for monitoring of the position of a single Abrikosov vortex 
in thermal depinning and Lorentz force depinning experiments. The recent experiments on 
a spontaneous fluxon formation in annular Josephson junctions demonstrate the necessity to 
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take into account the Abrikosov vortices trapping during the thermal quench as a compet- 



itive effect 



251 ] . In addition, a new detection principle based on the interaction of a single 



gamma-photon with trapped Abrikosov vortex is proposed for the development of a gamma 
-ray solid state detector with high intrinsic detection efficiency in the energy range up to 
100 keV |26|. 

In this paper we present a theoretical study of the Fiske steps for a "small" Josephson 
junctions in the presence of randomly distributed Abrikosov vortices. Our approach is 
based on the extension of a well known Kulik analysis for low-dissipative uniform Josephson 
junctions to the Josephson junctions with inhomogeneities. We obtain a peculiar regime 
where the Fiske steps amplitude shows a weak non-monotonic decrease with the vortex 
density n^. We experimentally trapped different number of Abrikosov vortices and for each 
vortex density ha we measured the dependencies of the critical current J c and the Fiske 
steps amplitude Ip on and on externally applied parallel to the junction plane magnetic 
field B\\. A good agreement between our theory and experimental results was found. 

The paper is organized as follows. In section [TT] the theoretical model of a Josephson 
junction with randomly distributed misaligned Abrikosov vortices is presented. In section 
IIHI by making use of a generic approach elaborated in the Ref . 27j , we calculate the Fiske 
resonances in inhomogeneous Josephson junctions, and in the section HVl the dependence of 
the Fiske steps amplitude Ip on the density of Abrikosov vortices is analyzed. Section 
IVl is dedicated to the experimental details, namely, the experimental setup, the sample 
description and the procedure to measure the Josephson critical current I c and the amplitude 
of first Fiske step Ip. In section IVTI we present the experimental results and the comparison 
with the theory. Section IVHI provides conclusions. 



II. MODEL OF A JOSEPHSON JUNCTION WITH RANDOMLY DIS- 
TRIBUTED MISALIGNED ABRIKOSOV VORTICES 

We consider a small, i.e. W < Aj, Josephson junction in the presence of randomly 
distributed pinned misaligned Abrikosov vortices (here, W is the size of a Josephson junction, 
and A j is the Josephson penetration depth). A magnetic field of a misaligned Abrikosov 
vortex enters and leaves the superconducting electrodes of a junction in different points. 
The distance between these points determines the misalignment length 5 (see Fig. 1). 
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FIG. 1: Schematic cross-section of a Josephson tunnel junction with a misaligned Abrikosov vortex. 
S is the misalignment length. A^ and Xl 2 are the London penetration depths of the superconduc- 
tors forming the junction. The orientation of the two magnetic field, B± and Bu, is provided. 



Under these conditions the Abrikosov vortex contribution <pv{p) t° the total Josephson 
phase difference ip(p,t) randomly depends on the position p in the plane of the junction. 



The particular 



in Refs. 



19 



2Q|. 



orm of a single Abrikosov vortex contribution has been obtained previously 



The dynamics of a Josephson junction with randomly distributed Abrikosov vortices is 
described by an inhomogeneous sine-Gordon equation jlj, [28 



d 2 <Pi(p,t) 1 d 2 <pi(p,t) <y d(p x (p,t) 



dp 



7& 



dt 2 



Of 



=A J 2 sin[yjy(p) +koX + u Q t + (pi(p,t)], 



(2) 



where 



k 



2tt j B||(A l1 + A L2 ) 



2eV 



$0 h 

Here V is the external dc voltage, cq is the velocity of Swihart waves, B\\ is an external 
magnetic field applied parallel to the y-axis, \l x and Al 2 are the London penetration depths 
of the superconductors forming the junction, the parameter 7 determines the decay of the 
electromagnetic waves in the junction and depends on the quasiparticle resistance of the 
junction in a subgap region, $ is the flux quantum . 
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III. FISKE RESONANCES IN INHOMOGENEOUS JOSEPHSON TUNNEL 
JUNCTIONS: THEORY 



Our theoretical analysis is based on the extension of t 
small low-dissipative Josephson junctions with 7 << u 
in the following form 



re Kulik approach developed for 
I27] ]. Thus, we present (pi(p,t) 



1. 



^(p, t) = J2 HeKe-*-*-*" *}, (3) 

n 

where the wave vectors of electromagnetic waves that can be excited in the junction, i.e. 
k n = nn/W, are determined by the junction's geometry. We note that in the particular 
case when the parameter 7 is not so small, Fiske resonances overlap each other, and a single 

n 

Eck peak is formed in the I — V characteristic of a Josephson junction Such feature 
has been theoretically analyzed also for Josephson junctions with randomly distributed 
inhomogeneities 



As the parameter 7 << uq the Fiske resonances are well separated, and for the analysis 
of n-th Fiske step we can choose a single wave vector k n in the Eq. (jSJ). Substituting Eq.flSJ) 
in Eq.(J2} and using the condition hco = c k n , i.e. ([T]), we obtain the nonlinear equation for 
the amplitude of Swihart electromagnetic waves as 

c o [J 2 (^) + Jffy]Re[ [ e i<pv(p)+ik *-ik n ^ (4) 



n A 2 j7 cV~ uv 2 ' ' ~ iV 2 /J_ " 7 S 
where Jo and J\ are the Bessel functions, S is the junction area. 
Similarly we obtain the expression for the amplitude of n-th Fiske step 

4 W) = JeJo^JiCYW/^e^^*-*^ . (5) 
where j c is the Josephson critical current density in the absence of vortices. 



IV. CRITICAL CURRENT AND FISKE RESONANCES: THE DEPENDENCE 
ON THE ABRIKOSOV VORTICES DENSITY 

Since the contribution of Abrikosov vortices tpv{p) to the Josephson phase difference is a 
random function of the coordinate p, we obtain the averaged quantities only. First we define 
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the correlation area for Abrikosov vortices as 



where the sign < .... > means the averaged value over the random positions of Abrikosov 
vortices. The correlation area <r(k) depends on the density of misaligned Abrikosov vortices 
ha, the misalignment length S, and the wave vector k that, in turn, is determined by an 
externally applied magnetic field parallel to the junction plane B». For the particular case 
of a circular junction the cr(k) has been found in 1 171 120| . and in the absence of B\\ 



a ~ Sexp(-nn A 5 2 \n(2L/5)) , (7) 

where L is the radius of the circular junction. It has been obtained in [if], [it} Q] that the 
critical current of the Josephson junction with randomly distributed misaligned Abrikosov 
vortices is determined by the correlation area cr{k) as 



J' 



V^kjS. (8) 



Thus, in the absence of B\\ as the critical current reaches the maximum, we obtain a strong 
decrease of the critical current I c with the Abrikosov vortices density ua'- 

laT = —' n 2L< (9) 
where J c o is the critical current in the absence of Abrikosov vortices. 

Next, we turn to the analysis of Eqs. fl4]) and (JSJ) determining the amplitudes of Fiske 
resonances. First, we notice that even in the absence of k (or B») the amplitudes of Fiske 
resonances 1^ are not zero. However, the magnitudes of Fiske steps can be increased by 
tuning of B» . The maximum values of Ip reach as the condition k n = ko is satisfied. In the 
following we will assume that the correlation area a is small, i.e. cr < 5. Using this assump- 
tion we express the maximum amplitudes of Fiske resonances 1^ through a single parameter 

er(0) = a. Indeed, in the limit as a is small or more precisely a < S(^p-f = a% the 

c o 

amplitudes a n are small, and expanding the Bessel function over a small argument a n we 
obtain 

r (i) _ jcP t [S 2( A , x-l 



If =jfc = jc<( AX pUor\ cr^a* (10) 
In this regime Ip is proportional to a and displays a strong exponential decrease with ua- 



In the opposite limit a > a*, the amplitude of electromagnetic waves becomes large but 
the oscillations of Bessel functions are strongly suppressed. Therefore, the amplitude of ijp 
is still small. The averaged value of Ip can be also expressed through a single parameter a 

as 

4 1 ' = (f ) 1/3 ^(^?) 1/6 )e-i [ - ]1/3 , o » a* (11) 

In this limit, the amplitude of Fiske resonance weakly depends on the density of Abrikosov 
vortices, displaying a small maximum on o ~ a*. Notice here, that this regime can be 
easily realized for low-dissipative junctions as 7 << uq. 



V. JOSEPHSON JUNCTIONS WITH MISALIGNED ABRIKOSOV VORTICES: 
EXPERIMENTAL SETUP AND MEASUREMENTS 

The experiments were carried out on the Nb — AlO x — Nb Josephson tunnel junction of 
square geometry and dimensions of 50 x 50 /im 2 . Details of the fabrication process are 



reported in 



29 



30] 



The junction was mounted on the sample holder inside a closed copper box in order to 
prevent the influence of the electromagnetic noise. The box was immersed in a liquid helium 
bath, so all measurements were done at a stable temperature T = 4.2K. In order to heat 
the junction above the superconducting transition temperature T c , a SMD 100 Ohm resistor 
was used as a heater mounted in contact with the sample holder. A special coil system 
containing two Helmholtz copper coils was used in order to generate both the magnetic fields 
perpendicular to the junction plane B± and the magnetic field parallel to the junction plane 
B\\ . The shielding of the sample from the Earth magnetic field and from the electromagnetic 
noise was reached by using a three level /z-metal shielding system and a one level /U-shield. 

The accurate measurements both the Josephson critical current I c , and the amplitude 
of the first Fiske step Ijp were achieved by a careful shielding and by inserting cold low 
pass filter, mounted just near the sample. First, we measured the I — V curves of the 
junction without trapped Abrikosov vortices. We tuned the amplitudes of Fiske resonances 
by an externally applied magnetic field JB11. Fig.2 shows the typical I — V curves where 
the amplitude of the first Fiske step was maximized (open circles) and suppressed (closed 
circles). The first I — V curve was used in order to obtain the voltage position of the Fiske 
step Vi„ while the second one was utilized for experimental estimation of the quasiparticle 
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FIG. 2: I — V curves recorded at B\\ = 2.2G when the first Fiske step is maximized (open circles) 
and at Bn = \2.1G when the Fiske step is suppressed (closed circles, ). The method allowing to 
extract both the amplitude of the first Fiske step, Ip \ and its voltage position, V\ is shown. The 
voltage position is obtained by the dashed line which is parallel to the current axis and intersects 
the Fiske step branch at I = l£ /2. The amplitude of the first Fiske step is 1^ = lj? — I q where 
I q is the quasiparticle current at voltage V\. V Ic is the threshold voltage to define the Josephson 
critical current. 

current in the sub-gap region. 

The method to determine the voltage position V\ of the first Fiske step is shown in Fig. 2. 
We note that the voltage position can also be estimated as the voltage at current If (see 
Fig.2) but the difference between this method and the one represented in Fig.2 is negligible. 

The experimental technique to measure the Josephson critical current I c was as follows. 
The current / was increased linearly in time from zero value, and the junction voltage V was 
simultaneously recorded. When the current reached the Josephson critical current value, a 
non zero junction voltage appeared. Under a condition that the junction voltage equals to 
the threshold voltage V Ic (criterion of the Josephson critical current, see Fig.2), the current 
value was recorded and the current returned back to zero and then the next measuring cycle 
started again. The current sweep frequency was 10 Hz and the value of V Ic was of 10 jiV . 

A cycle to measure the amplitude of the first Fiske step was similar to the Josephson 
critical current measurements but the starting value of the current was chosen to satisfy the 
condition that the junction voltage was greater then zero and less then voltage position of 
the first Fiske step. From this starting point the current was increased, the first Fiske step 



9 



branch was recorded and the maximum current at the first Fiske step branch 1^ was defined 
as the current at which the junction voltage became greater then V\ + V Ic (see Fig.2). The 
amplitude of the first Fiske step was determined by substracting the quasiparticle current at 
Fiske step voltage I q from the switching current at Fiske branch 1^ . Both Josephson critical 
current and amplitude of the first Fiske step for a given value of parallel magnetic field Bn 
were obtained. The I C (B\\) and Ip\B\\) dependencies were measured for both polarities of 
the bias current. 

The experimental cycle for the evaluation of the influence of Abrikosov vortices on the 
Fiske step amplitude consisted of the following steps: 

1) The sample was heated to T > T c ; 

2) A magnetic field B± of fixed value was applied; 

3) The junction was cooled in the field B± up to T = 4.2 K, i.e. immersed in the liquid 
helium at normal pressure; 

4) The magnetic field B± was turned off; 

5) The amplitude of the first Fiske step was maximized by the parallel magnetic field B\\ 
and the I — V curve with the first Fiske step branch was measured in order to obtain the 
voltage position of the first Fiske step; 

6) The I C (B\\) curve was measured; 

7) The I^{B\\) curve was measured; 

8) Both the Josephson critical current and the amplitude of the first Fiske step, were sup- 
pressed by a suitable parallel magnetic field B\\ and the I — V curve was recorded in order 
to measure the quasiparticle current at V — V\. 

After the field cooling process the Abrikosov vortex density can be estimated as 

n A ~ (12) 

VI. EXPERIMENTAL RESULTS AND DISCUSSION 
A. Junction Parameters 

The measurements of the I C {B\\) dependence without trapped Abrikosov vortices were 
performed at T = 4.2 K (Fig. 3, open circles). The Josephson penetration depth \j was 
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FIG. 3: The experimental I C (B«) curve (open circles) and theoretical fit by Eq.(10) (solid line). 
The trapped Abrikosov vortices are absent. 

estimated to be about 54 fim. The experimental I C {B\\) curve is in a good agreement with 
the theoretical one which was calculated for a "small" square junction by the Fraunhofer 
formula [l| 



1,(5,1) = ]c W 2 



wBuWd 



(13) 



*0 

where j c is the Josephson critical current density, W is the dimension of the junction, 
d = A/ /1 + Xl 2 + t and Al x and Xl 2 are the effective London penetration depths of the 
superconductors forming the junction, and t is the barrier thickness. The theoretical curve 
calculated by Eq.( TT3~|) is reported in Fig.3 by solid line. This good agreement between 
experimental and theoretical results confirms the uniform distribution of the Josephson 
current density over the junction area. 

B. Magnetic field dependence of the first Fiske Step without Abrikosov vortices 

In the absence of trapped Abrikosov vortices the value of voltage position of the first 
Fiske step V\ was determined from the I — V curve at T = 4.2 K when the Fiske step 
amplitude was maximized by magnetic field B» (Fig. 2, open circles). Using this value of V\ 
and Eq.(l), we calculated the resonant frequency loq and Swihart velocity Cn. The values of 
V\, ojq and Co are given in Table 1. 

The ijp (fi,|) curve measured after cooling the junction in a zero magnetic field B± (open 
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TABLE I: Experimental values of the voltage position of the first Fiske step V±, the Swihart velocity, 
cq, and the resonant frequency, u>q, obtained from Eq.(l). 



Vi (jjV) 


c (m/s) 


u (s x ) 


226.3 


1.1 x 10 7 


6.9 x 10 11 




Parallel Magnetic Field (G) 



FIG. 4: Experimental iP(Bh) curve for the junction without trapped Abrikosov (open circles). 
Solid line represents the theoretical dependence calculated by Eq.(|3J) and Eq.© for 7 = 1.1 x 
10 10 s _1 . Dashed line is the result of calculation without J%(%) term at 7 = 4.3 x 10 9 s _1 according 
to the original theoretical analysis of [271 ] . 

circles), is shown in Fig. 4. The theoretical dependence for the first Fiske step amplitude 
on field B» was numerically computed from Eq. (j4j) and Eq. (jSj) in the absence of Abrikosov 
vortices. The value of the parameter 7 = 1.1 x 10 10 s -1 provides a best fitting of ex- 
perimental data. Similar value of 7 can be extracted from the I-V curves where the Fiske 
resonances were suppressed. The agreement with the theoretical analysis is excellent even 



for the second and third lobes. In the Ref. 



27J a similar theoretical analysis of the Fiske 



resonances of a junction without Abrikosov vortices has been carried out. However, in the 
equation determining the amplitude a of electromagnetic waves the term of J\ was omitted. 
Taking into account this term one can get a better agreement with experimental data. 
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FIG. 5: I c (-B||) curves measured after cooling of the junction in perpendicular magnetic field B± 
of various values. 

C. Josephson critical current of a junction with trapped Abrikosov vortices 

Fig. 5 shows I C {.B\\) curves measured after cooling of the sample in different perpendicular 
magnetic fields B±. We obtained from each I C (B\\) curve the maximum values of I c for both 
polarities, I+ max and I~ max , and then we selected the maximum value between I+ max and 
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FIG. 6: /> j (S||) curves measured after the junction was cooled in perpendicular magnetic field 
H per of various values. 

j-max d eno ting it as I™ ax . Similar to H]' we found the deviation of the I C (B\\) curve 
from a "vortex-free" Fraunhofer-like curve (Fig. 3, open circles) and we observed a strong 
suppression of I™ ax . For each value of B± the vortex density ua was calculated by using 
Eq. (|T2l) . The dependence of the maximum Josephson critical current l™ ax on the density 
of trapped Abrikosov vortices is shown in Fig. 7 by open squares. It is well known that 
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during the field cooling process, misaligned Abrikosov vortices are randomly trapped over 
the junction area (see Fig.l). Eq.flU]) describes the dependence of I™ ax both on the vortex 
density and on the misalignment length 8 in case of a circular junction. Taking L to 
be the radius of a circular junction with an area equal 50 x 50 fim 2 , we used Eq.([H]) and 
Eq. ( II 2p in order to extract the value of 5 for each value of vortex density. The average value 
of misalignment parameter over the all vortex densities is 5 = 0.8 ± 0.1 fim. 



D. First Fiske step in the presence of trapped Abrikosov vortices 

We observed that the presence of trapped Abrikosov vortices leads to a great deviation 
of the I { f\B\\) curves from the "vortex- free" one shown in Fig. 4 by open circles. These 
curves are shown in Fig. 6. The presence of trapped vortices did not lead to the variation of 
the voltage positions of Fiske resonance. This result indicates that the trapped Abrikosov 
vortices did not change the geometrical conditions for appearance of the first Fiske step and 
only influenced on the magnitude of Iff 1 by additional spatial variation of the Josephson 
phase difference ip(p,t). 

Next, we notice that a non-zero amplitude of the first Fiske step at Bn = appeared as 
Abrikosov vortices were trapped. However, by tuning the parallel magnetic field we could 
substantially increase the amplitude of the Fiske step. Using the same procedure which 
was used to obtain the I™ ax from the I C (B\\) curves (see subsection IVI CI ), we found the 
maximum amplitude of the first Fiske step j^' max f or each value of B±. The experimental 
dependence of j^' max on the Abrikosov vortex density is shown in Fig. 7 (open circles). 
In contrast to the J™ ax (n^) dependence, the weak nonmonotonic decrease of the magnitude 
of I { p> max with tia was observed. Such behavior of the I^ 1,max {riA) dependence indicates 
that the condition of a > a* should be applied for our sample. Indeed, this conclusion is 
confirmed by estimating the values of both a and cr* for our junction. 

In order to estimate the critical value of cr*, we have to obtain the dissipation parameter 
7. This parameter determines the decay of the electromagnetic waves in the junction, and it 
can be obtained as 7 ~ {R q C)~ where C is the capacitance of a junction per unit area, R q is 
the quasiparticle tunneling junction resistance at voltage V = V\ per unit area [l|. T rapped 
Abrikosov vortices increase the quasiparticle current by virtue of their normal cores 
and, the value of 7 has to also increase. Therefore the value of 7 = 1.1 x 10 10 obtained 
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FIG. 7: Open squares: Experimental dependence of the maximum Josephson critical current /™ ax 
on the density of trapped Abrikosov vortices n; Open circles: Experimental dependence of the 
maximum amplitude of the first Fiske step jW> maa: on the density of trapped Abrikosov vortices 
n. Solid line is the theoretical Ip^' max (ua) curve, calculated by Eqs. pTj) and (|7|). The values of 
7 = 2.2 x 10 10 s -1 and a mean misalignment parameter over all vortex densities 5 = 0.8fim were 
used; Full circles represent the theoretical Ic 1,max (tia) dependence obtained by Eqs. lfTT]) and ([8]). 

for the junction in the absence of Abrikosov vortices, was not convenient for the case of 
trapped Abrikosov vortices. We used 7 as a fitting parameter in order to reach the good 
agreement between the experimental value of 1^ recorded at Bj_ — 1.1 G and the theoretical 
value calculated by using Eq. fll]) and Eq. ([5]) for the same magnetic field. We obtain 
7 = 2.2 x 10 10 s -1 from this fitting procedure and the parameter cr* = 3.2 x 10 2 [im 2 was 
estimated. It is important to note that the quality factor Q = uq/^ oi the Fiske resonance 
is equal to 31, such that a crucial assumption of low dissipative junctions is valid for our 
junctions. 

Using Eq. ([7]) we calculated the correlation area a for each value of density tia, and we 

16 



obtained that for all values of B±, the correlation area a > cr*, and therefore, the Eq. (ITT|) 
is valid. The theoretical dependence of Ip on ua is shown in Fig. 7 by solid line. 

The experimental dependence Ip^' max (ua) agrees qualitatively with the behavior of theo- 
retical curve but several experimental points disagree with our theoretical predictions. One 
of possible reason of this discrepancy is the variation of the misalignment length 6 for different 
magnetic fields B±. In order to check this hypothesis we extracted a from the experimental 
values of I™ ax by using Eq. (JSJ) for different values of ua- Substituting the obtained values 
of cr into Eq.flTTJ we calculated j^' max f or each experimental values of ua- The result of 
such procedure is shown in Fig. 7 by closed circles. By making use of this procedure the 
agreement between the theoretical dependence and experiment becomes better. So we can 
conclude that the variation of 5 strongly influences the maximum amplitude of the first Fiske 
step. 



However, we have to notice that our 



only. As has been shown in Refs. 20 



theoretical analysis is devoted to averaged quantity 



211 ] . in inhomogeneous Josephson junctions strong 



mesoscopic fluctuations of various physical quantities occur. These fluctuations manifest 
themselves in the form of random oscillations on the dependence of Ic{ua) (or I C (B\\\). 
These oscillations can be especially large for the dependence p^' max on nA d ue to a strong 
nonlinearity of Eqs. @ and (JSJ). 

The pronounced deviation of the experimental point at tia = 0.15 fim~ 2 (Bj_ = 3.2 G) 
and ua = 0.26 ^m~ 2 (B± = 5.3 G) from the theoretical prediction (see Fig. 7) can be also 
attributed to the trapping of so called "monopole" Abrikosov vortex. The "monopole" vortex 
penetrates only a one superconducting electrode of the tunnel junction and its magnetic field 
spreads over the total area of a "small" Josephson junction. Such Abrikosov vortex trapped 
near the junction center, alters the I C {B\\) dependence from the Fraunhofer pattern to the 
one with two primary lobes and with suppressed Josephson critical current at B» = 0. 
In addition, the amplitude of the first Fiske step is maximized. Even if a large number of 
Abrikosov vortices is trap ped , the effect of the monopole vortex, pinned in the central part of 
junction area, dominates [32|. Because of a statistical character of the trapping process, the 
cooling in the fields B± = 3.2 G and B± = 5.3 G can provide the trapping of "monopole" 
Abrikosov vortex near the junction center, which, in turn, maximize the first Fiske step, 
similar to the magnetic flux trapping effect in case of a Josephson junction with annular 
geometry 33[. The shapes of the I C (B\\) dependencies shown in Fig. 5c and 5e (especially 
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the curve e), permit to sustain our explanation of the deviation of the experimental data 
from theoretical prediction presented in Fig. 7. 



VII. CONCLUSIONS 

We have developed a theoretical and experimental study of the Fiske step resonances 
in low-dissipative Josephson junctions with randomly distributed misaligned Abrikosov vor- 
tices. This analysis is an extension of the approach elaborated in Ref. [27| to inhomogeneous 
Josephson junctions. The local inhomogeneous magnetic field of misaligned Abrikosov vor- 
tices leads to an additional contribution to the Josephson phase difference. The influence of 
randomly distributed misaligned Abrikosov vortices was described in terms of the correlation 
area a (see, Eq. (JS}) which depends both on the Abrikosov vortex density ua and on the 
misalignment length of vortices 5. We obtained a peculiar regime as a > a* = S( Xj1 ^ ) 2 
characterized by a weak non-monotonic dependence of the maximum Fiske step amplitude 
on ua- In this regime the amplitude of excited electromagnetic waves in the junction a n > 1. 
This behaviour is in a sharp contrast with a strong decrease of the maximum critical current 
I c with riA- The critical parameter a* depends on the dissipation 7 in the junction. For 
low-dissipative junctions the parameter a* << S. 

We carried out the experiments on small low-dissipative Nb based Josephson tunnel 
junction with artificially introduced Abrikosov vortices. We measured both the critical 
current I c and the Fiske resonant steps in the I — V curves. We observed that the critical 
current I c strongly decreases with the density of Abrikosov vortices in a complete agreement 
with a theory, Eq. (JHJ). 

In the absence of trapped Abrikosov vortices the experimental dependence of the ampli- 
tude of first Fiske resonance on B\\ is in a good agreement with the theoretical description 
by Eqs. @ and flH]), where the term J\ is taken into account. As the misaligned Abrikosov 
vortices were introduced in the junction, Fiske resonances were observed even in the absence 
of B\\. However, an application of a small magnetic field B\\ allowed to increase the ampli- 
tude of Fiske resonances. The voltage position of the first Fiske resonance determined by a 
geometry of the junction, did not vary in the presence of Abrikosov vortices. The measured 
dependence of the maximum amplitude of first Fiske step j^' max on nA showed the qual- 
itative agreement with our theoretical analysis (Eq. QTjl ) i.e. with a weak non-monotonic 
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dependence (see Fig. 7). Therefore, the regime a > a* for our low-dissipative junctions was 
realized. 

Finally, we notice that Eq. ( TTTj) describes the averaged value of the amplitude of Fiske 



resonances. However, as was indicated in Refs. |20l . |21| . strong mesoscopic fluctuations 
of various physical values can be observed in inhomogeneous Josephson junctions. These 
mesoscopic ffuctuiations are especially strong for the amplitude of Fiske resonances due to 
nonlinear character of Eqs. (j4]) and fl5]). Therefore, observed deviations of some experimental 
point from the theoretical curve described by Eq. lflT]) . can be attributed to such mesoscopic 
fluctuations. This discrepancy could be also due to a possible trapping of a single monopole 
Abrikosov vortex just near the junction center. 
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